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ABSTRACT 

We consider the radial buoyancy of vertical magnetic field lines in radiation 
and gas pressure dominated accretion disks. We find that in addition to radial 
drift driven by turbulent diffusion and biased by the global field geometry, 
there are buoyancy effects which tend to move magnetic flux outward. In gas 
pressure dominated disks the poloidal magnetic field will move outward at a 
rate comparable to its inward advection. On the other hand, in a radiation 
pressure dominated disk the poloidal magnetic field will usually move outward 
faster than it is advected inward. This implies that the fields in disks in active 
galactic nuclei are generated at small radii by an internal disk dynamo. This 
conclusion can be avoided if the external field imposes a supersonic Alfven speed 
within the disk without giving rise to interchange instabilities. In any case we 
note that variations in the mass transfer rate will lead directly to a modulation 
of the nonthermal emission from the disk system. 
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1. Introduction 

It is well known that the emission- line spectra of active galactic nuclei (hereafter 
AGNs) are strongly variable. Typical time scales for this variation range continuously from 
the time scales appropriate for light travel across the inner regions of the accretion disk 
around a black hole, i.e. 

t « GM/c 3 , (1) 

(hours for supermassive black holes) up to years. This variability remains one of the central 
issues in any theoretical model of the central engine in an AGN. 

A popular model for the source of the nonthermal radiation from an active galactic 
nucleus (AGN) consists of a supermassive black hole surrounded by a magnetized accretion 
disk. In its simplest form this model is taken to be axisymmetric and stationary. The 
nonthermal radiation is ascribed to the Blandford-Znajek process (Blandford and Znajek 
1977), which emits nonthermal radiation from an extended magnetic field embedded in 
the black hole event horizon, the surrounding plasma, and the accretion disk. This model 
was reformulated and extended by Thorne and Macdonald (1982), Macdonald and Thorne 
(1982), and Thorne et al. (1986) in the '3+l'-spacetime formalism. The magnetic field in 
the disk will also give rise to a magnetically driven wind, and possibly a jet, as matter 
spirals out along the poloidal field lines (cf. [Blandford fc Payne 1982j , [Najita fc Shu 199~4| , 



[Shu, Najita, Ruden, fc Lizano 1994 , and |Shu, Najita, Ostriker, Ruden, Wilken, fc Lizano 



|1994|) . In addition, the ejection of toroidal flux from the disk may lead to winds and jets 



e.g. |Contopoulos 1995]) 



Starting from the axisymmetric, stationary model for emission due to the Blandford- 
Znajek process Park and Vishniac (1989a, 19896) explored a time- dependent model which 
included the effects of variations in the mass accretion rate. The main point was to add 
the secular effects of mass accretion to the original axisymmetric, stationary model of 
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the nonthermal radiation. Such variations are expected because AGN accretion disks 
are expected to be radiation pressure and electron scattering dominated, and this can 



give rise to thermal and viscous instabilities on a broad range of time scales ([Lightman 
|fc Eardley 1974j) . In these two papers we investigated the axisymmetric, nonstationary 
electrodynamics of a black hole and its accretion disk, respectively. 

In the latter paper we found that the electrodynamic power output from the accretion 
disk can be variable on time scales associated with secular disk instabilities. The point was 
that the local fluctuations in fluid velocities in the accretion disk will cause fluctuations in 
the nonthermal component of the radiation by the Blandford-Znajek process. The time 
scales for these fluctuations, therefore, will reflect the range of orbital periods in the inner 
annulus of the disk. 

In this paper we discuss another effect which may explain the fluctuations in the 
nonthermal emission, the relative motion of the infalling matter and the entrained poloidal 
field. The traditional view of the magnetic field in astrophysical objects, like accretion 
disks, is that the matter and magnetic field lines are tightly anchored to one another. 
Consequently accretion disks should continuously accrete poloidal field from the surrounding 
medium. Since the back reaction from the field on the matter flow will also increase without 
limit this leads to a situation in which the magnetic field completely overwhelms the gas 
pressure and drastically alters the properties of the accretion disk. 

On the other hand, if one supposes that the magnetic field is not tightly coupled to the 
underlying fluid, then one can invoke the concept of turbulent diffusion for the magnetic 
field. This leads to the prediction that large scale magnetic fields tend to move outward in 
accretion disks and are eventually ejected from their outer edges ( |Van Ballegooijen 1989| ). 
Although it is difficult to reconcile turbulent diffusion with the view that the magnetic field 
fills the entire fluid volume, it is consistent with the idea that the magnetic field in stars 
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and accretion disks naturally forms intense, partially evacuated, fibrils that are capable of 
moving relative to the bulk of the fluid, i.e. flux tubes. In a recent paper Vishniac (1995a) 
examined the underlying physics behind the formation of such flux tubes and compared the 
predictions of a particular model with numerical simulations. A subsequent paper (Vishniac 
1995b) predicted that the flux tubes in ionized accretion disks will be largely empty and 
discussed the consequences of this model for the vertical buoyancy of the toroidal field of 
an accretion disk. Here we extend this work to the radial motion of a poloidal field and 
show that in a radiation pressure dominated disk it leads to a strong outward motion of 
the field lines. This in turn leads to a modulation of the nonthermal power emitted by the 
magnetosphere as changes in the state of the disk raise and lower the rate at which the 
poloidal field is dispersed outward. 

In §11 of this paper we discuss the radial buoyancy of vertical magnetic field lines 
embedded in an accretion disk. In §111 we discuss some of the implications of this work for 
nonthermal radiation from AGN. 



2. The Radial Buoyancy of B z 

In order to understand the transport of flux within an accretion disk we need to have 
some model for the turbulent motions with the disk. Our work here is based on the idea 
that in an accretion disk these motions are driven by the magnetic field with the result 
that the angular momentum transport takes place within the disk and is mediated both by 
hydrodynamic transport and shearing of the local magnetic field. A magnetized accretion 
disk will be subject to a violent shearing instability which will act to transport angular 
momentum outward ([Velikhov 1959| , |Chandrasekhar 1960| , [Balbus fc Hawley 1991|) . The 



corresponding dimensionless viscosity a will be roughly the ratio of magnetic to thermal 
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B " ( Va \ 

a ~ TTi ~ — ' ( 2 ) 



energies or 

a ~ 

8ttP V c s 

where B and P are the characteristic values of the magnetic field and pressure, respectively, 
in the disk, and Va and c s are the corresponding values of the Alfven speed and sound speed. 
The resulting turbulence is anisotropic, but only by factors of order unity. This turbulence 
will be characterized by an eddy size of order Va/Q and a typical turbulent velocity of order 
Va- Larger scale modes, as well as the Parker instability, will be largely suppressed by these 
eddies flVishniac fc Diamond 1992 ). Nevertheless, there will be a residual mean vertical 



buoyant velocity of order ac s . A similar result applies when the magnetic field is considered 
to consist of individual flux tubes instead of some diffuse, space-filling mean field QVishniac" 



1995b|) . For a thin disk the sound speed c s ~ Htt, where H is the disk thickness and f2 is 



the rotational frequency. Consequently, the dominant turbulent eddy size will be 

L T ~ a 1/2 H. (3) 



The magnetic field in a turbulent medium will tend to separate into regions of high 
and low magnetic field density. Most of the magnetic energy will be contained in flux tubes 
whose internal density is kept at a level well below that of the ambient medium through 
turbulent pumping QVishniac 1995a|) . These flux tubes will interact with the surrounding 
fluid through the turbulent drag of the fluid as it moves by. Their mutual interactions 
guarantee that they will be broadly distributed in the fluid. The force per unit length on a 
flux tube is 

Pr^M) = Cd P^ T ~ ^\(Vt ~ V t )r t + n c (nr])^- - (p - Pi )^7rr t 2 ), (4) 

where the subscript t denotes quantities evaluated within the flux tubes, R c is the radius of 
curvature of the local magnetic field, Vr refers to the local fluid velocity, is the coefficient 
of turbulent drag (of order unity), g is the local gravitational acceleration, and h c is a unit 
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vector pointing in the direction of the curvature of the flux tube. We will subsequently 
use \Vt\ to refer to the root mean square turbulent velocity in the fluid. Only the velocity 
components perpendicular to the local magnetic field enter into equation If the flux 
tubes are largely empty, as they are in highly ionized disks, then the right hand side of this 
equation must sum to zero. Under most circumstances the buoyant force term will be much 
smaller than the other two terms ( |Vishniac 1995 b| ) . The properties of the flux tubes are 



largely independent of the resistivity of the gas provided that 

«(i§)(ST"">(4)"' 

where Air is the turbulent Mach number, is the wave number of the large scale eddies, 
1] is the resistivity, Eb and Et are, respectively, the spatially averaged magnetic and kinetic 
energy densities, f c , and w are constants of order unity, 7 is the adiabatic index of the fluid, 
and n is the turbulent power spectrum index (equal to 2/3 for Kolmogorov turbulence). 
The viscosity can be ignored if 

These conditions define the ideal fluid regime, within which the flux tubes are largely 
evacuated, except for a surface layer whose width is proportional to the square root of 
the resistivity. In accretion disks the kinetic and magnetic energy densities are roughly 
comparable due to the Balbus-Hawley instability. In this case, equations (j5|) and (H) can be 
rewritten as ( |Vishniac 1995b| ) 



and 



., :c«h\ ( p 

120 <~ a' 



v \P 

u ' \ 1 gas , 



where P is the total pressure, and P gas is the gas pressure. For temperatures and Mach 
numbers appropriate for ionized disks, both of these conditions are usually satisfied. In this 
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regime the typical flux tube radius is ( |Vishniac 1995a] , |Vishniac 1995b| ) 



r t ~ o?' 2 H^-. (9) 

"gas 

The factor of P/P gas enters here because turbulent pumping is ineffective for photons 
in accretion disks. Consequently, the magnetic pressure in a flux tube is limited by the 
ambient gas pressure instead of the total pressure. 

In order to calculate the effects of radial buoyancy, we need to reconsider the third 
term in equation (^). The disk properties will vary radially, and the flux tube's length and 
radius will change as it moves. We need to calculate the change in energy associated with 
moving a flux tube radially. If we consider an empty flux tube of radius r t and length L 
embedded in an accretion disk with an ambient pressure P a length L then it will displace 
a volume of fluid equal to its own volume. Removing (adding) such a flux tube from an 
accretion disk will cause the surrounding gas to expand (contract) and lead to a drop (rise) 
in the local energy due to pressure work. The energy involved is 



Utube — A 



/oo 1 1 
(-pz 2 n 2 + -P)dz 
-oo I 7—1 



(10) 



7 

If we assume that the gas expands adiabatically to fill the volume of the flux tube as it 
moves, then we can use the virial theorem for a pressure supported thin accretion disk, i.e. 



/OO /"OO 
pz 2 tt 2 dz= / Pdz, (11) 
-oo J —oo 

to show that 

/oo AT 
P{z)*r\{z) a —dz ~ P(nr 2 )L, (12) 
-oo CtZ 

where we take the pressure and flux tube radius to be the appropriate values near the 
midplane of the disk. In other words, we obtain the usual expression for the energy loss 
due to adiabatic expansion. The contribution from the gravitational energy term, which is 
of comparable order, is canceled by an additional thermal energy term. This approach fails 
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for a flux tube in a radiation pressure dominated environment because only the gas pressure 
contributes to the energy expended by the gas expanding into the volume previously 
occupied by the flux tube. On the other hand, the gravitational potential term depends 
only on the matter density, and is not reduced by the tendency of the photons to fill up the 
otherwise empty flux tube. Consequently, U tu be is of order PrfL rather than P gas r 2 L. 

Moving a flux tube radially implies a change in U tu be and a corresponding radial 
buoyancy. The radial force on a flux tube is equal to the radial derivative of Utube- Since 
a given flux tube can split, or combine with other flux tubes, as it moves outward we need 
to constrain the radial variation of Utube by using the conservation of magnetic flux, i.e. 
B t r 2 oc Pgl^rf is a constant of motion. Consequently, the radial force per unit length due 
to buoyancy is 

F r » -P(nr 2 t )d r HPP^L). (13) 

Since the magnetic flux tube threads a turbulent disk, its actual length is 
~ H(H/L T ) ~ Ha- 1 ' 2 . This implies that 

F r « -P{-nr 2 t )d r l - \n{P 2 H 2 P; a \a' 1 ). (14) 

The sign of this radial derivative is not obvious, and may not even be the same for all radii 
and at all times. For a gas pressure dominated disk we have P = P gas and M w aPHQ -1 . 
Consequently, 

d r l - HP 2 H 2 p- a l) = d r l - \n(Mc s a~ 2 ) : (15) 

and the buoyant force will point outward in a stationary disk. On the other hand, near 
a thermal transition front, with d r T and d r M negative, magnetic flux lines can be pulled 
inward. 

For a radiation pressure dominated disk M « aPHQ-\ P gas oc EH^P 1 ' 4 , and 



-10- 



H oc M, so that 

oc T 7 M 3 Q. (16) 

P ry 

For a stationary disk this increases inward, so that once again we obtain an outward 
buoyant force on the magnetic field lines. 

If we assume that the curvature of a flux tube across the disk is negligible (more on 
that later) we obtain the mean radial motion due to buoyancy by equating the radial 
buoyant force given above with the turbulent drag due to a systematic outward motion of a 
flux tube. Averaging the turbulent drag on the flux tube over many eddies we find that 

(\V T - V t \(V T - V t )) « ~{\V T - V t \)V b ~ -\V T \V b , (17) 

where we have assumed that (Vt) = and (Vt) = V- Dropping constants of order unity we 
obtain 

— ~ p\V T \V r r t — . (18) 

This implies that 

V r - ^-a^. (19) 

1 gas 1 J L 

In other words, the radial buoyancy of the field lines gives rise to an outward drift which 
is of order the inward drift of the matter times P/P gas . This same procedure (balancing 
buoyancy with turbulent drag) was used in Vishniac (1995b) to obtain a vertical buoyant 
velocity of ~ ac s . Note that we have ignored the specific angular momentum of the gas 
entrained on the field lines. This is justified as long as the flux tubes are largely empty, as 
they are in ionized accretion disks. In a disk dominated by gas pressure, the inward moving 
matter will act to drag the field lines inward, and the actual motion of the field lines will 
be determined by the detailed balance between these two effects. However, in a radiation 
pressure dominated environment the field lines will usually move outward. This conclusion 
will be affected by viscous and thermal instabilities in AGN disks, which may temporarily 
produce radial gradients which will pull the magnetic field lines inward. 
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When the flux tubes are not empty they will exchange matter, and consequently 
angular momentum, with their environment at the eddy turn over rate, which is ~ Q. 
Neglecting other effects which will tend to equalize the specific angular momentum of 
material inside and outside flux tubes, we see that this exchange rate implies that the 
average specific angular momentum of the flux tubes can differ from the surrounding gas by 
as much as 

A(r 2 n tube ) ~ ^d r (r 2 n gas ). (20) 

Of course, there will be large fluctuations in this quantity, since the magnetic field will 
drive an instability by creating pockets of gas and field with specific angular momenta 
both larger and smaller than the local average. This systematic bias in the specific angular 
momentum of gas in the flux tubes will lead to a radial force of ~ —p t QV r . Comparing this 
to the buoyant force per unit volume and using equation ([IUp we have 

7 VY1 ~^ a - 21 

(P ~ Pt)ci/r p~ p t Pgas 

Since the radiation pressure differential across a flux tube is negligibly small ( |Vishniac 
1995 b|) the magnetic pressure in a tube must be balanced by a drop in the matter pressure. 



Therefore, 



(B 2 ) B\ AP gas (p - Pt ) P gas 

8vrP 8vrP P p P y 1 



From this we conclude that the ratio of the forces due to the specific angular momentum of 
material inside flux tubes to the buoyant forces will be less than pt/ p, which is less than one. 
Of course, this limit will be approached only when the flux tubes are barely distinguishable 
from the rest of the fluid, so that the volume average magnetic pressure is close to the 
magnetic pressure in the flux tubes. If we can neglect specific angular momentum, then the 
radial velocity due to buoyancy depends on the flux tube radius and internal magnetic field 
only in the combination B 2 r t , which is always of order V^L? for a magnetic field embedded 
in a turbulent medium and in equipartition with the turbulent energy ( IVishniac 1995a]) . 
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We conclude that the results presented here should apply to all accretion disks except for 
those where the resistivity is so large as to nearly obliterate all flux tube structure. 

Up to now we have assumed that the magnetic field outside the accretion disk crosses 
the disk with a negligible bending angle. This will not be the case if the external field is 
driving a magneto-centrifugal wind (e.g. piandford fc Payne 1982|) . In this case the motion 



of the field can determined by vertical mixing through the disk. We can see this effect for 
a weak external field by using the usual mean field theory, and ignoring the division of the 
magnetic field in the disk into flux tubes. Our argument is equivalent to the one given 
by Van Ballegooijen (1989). The only conceptual difference is that we have justified the 
concept of turbulent diffusion by appealing to a specific physical model for the microscopic 
structure of the magnetic field. Starting from the induction equation 

8 t B = V xvx B, (23) 

— * — * 

we can define a transient field b which is the response of the large scale magnetic field, B , 
to small scale turbulent motions. We have 

d t b = VxvxB . (24) 

Substituting this back into the right hand side of the induction equation we obtain 

9,5„ = Vx (vx (J v-VB dt)) -V x (v x (fv-VBdt)), (25) 

for incompressible turbulence. The first term in this equation is the source of the turbulent 
dynamo, if any, and the second term gives rise to diffusive effects. If we consider only the 
diffusive terms which affect the evolution of B z then we get 

8 t B z = -d r Uv 2 z r)d z B r ) + d r Uv 2 r r)d r B z ) + d r ((v z v r r)(d z B z - d r B r )) , (26) 



where r is the correlation time of the turbulence, in this case ~ Q 1 . The last term will be 
small, not only because the off diagonal terms of the averaged stress tensor will be small, 
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but because the extra radial derivative of B r will make it of order H/r compared to the 
first term. (Since the divergence of the mean magnetic field vanishes, the two parts of the 
last term are equal to each other.) This implies an outward drift of the vertical field lines 
given by 

V r = (v 2 z t) %^ ~ ac s tan 9, (27) 

where 9 is the bending angle of the magnetic field as it passes through the disk. Similar, and 
more precise, estimates can be found in paper by Lubow et al. (1994) and Reyes-Ruiz & 
Stepinski (1996). A different result can be obtained by assuming that the level of turbulent 
diffusivity in the disk is much less than that required to explain radial angular momentum 
transport. This is conceivable only if the Alfven velocity in the disk is supersonic, thereby 
suppressing the Balbus-Hawley mechanism. If the magnetic field is driving a wind, then 
9 will be positive and of order unity and this will dominate over the radial buoyancy of 
field lines for (P gas /P) > (H/r). On the other hand, if the field's geometry is not directly 
affected by the disk, and has a curvature scale comparable to r then 9 will be of order Hj r 
and the resulting radial drift velocity will be of the same order as the inward velocity of the 
matter. In this case radial buoyancy effects will dominate for radiation pressure dominated 
disks. Finally, in the event that the magnetic field lines bend inward as they cross the disk, 
this diffusive effect will move the vertical field lines inward. 

We note that this derivation neglects the tension due to the bending of the large scale 
field as it crosses the disk and the possibility that the angular momentum in the disk is 
transported out into a magnetically driven wind instead of radially within the disk. Both 
assumptions are justified for a weak external field. A precise estimate of when this will run 
into problems must necessarily await a consensus on the degree to which an external field 
can produce a torque. Successful wind and jet models based on an external poloidal field 
( IBlandford fc Payne 198^ , [Najita fc Shu 1994| , |Shu, Najita, Ruden, fc Lizano 1994 and 



[Shu, Najita, Ostriker, Ruden, Wilken, fc Lizano 1994Q give an equivalent 'a' based on an 
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external torque which can be as much as 

a ext ~ (28) 

but these models assume that large radial bending angles near the disk are possible. As 
mentioned above, this is reasonable only for an induced magnetic pressure in the disk which 
is at least as large as the ambient thermal pressure. At smaller values the external field will 
induce turbulence, and transport, within the disk. 

Without a definite model for the external torque, it is difficult to compare external 
and internal torques, but there are two points to consider with respect to the relationship 
between the internal turbulent viscosity and B ext . First, the formation of flux tubes implies 
that the same amount of flux is compressed into a smaller volume, giving rise to a higher 
internal rms B z than B ext . If the flux tube formation process is very efficient, as would 
expect in an ionized disk, then 

B Zi i n ternal ~ B]j^ xt (8irPg as ) (29) 

Moreover, the ai nt associated with an imposed vertical field is ~ (V^ 2 /c s ) rather than 
~ {Vaz/cs) 2 QVishniac fc Diamond 1993|) . Consequently, we have 

^fef ■ (30) 

Of course, this limit neglects any magnetic field generated within the disk itself, but such 
a field will only increase atint- We see that in spite of the potential advantage inherent 
in external torques due to their extended moment arm, reflected in the factor of r/H in 
equation (ESI), internal torques are likely to dominate for Va z < c s . 



Finally, we note that the second term in equation (^Bj) will tend to smooth out radial 
gradients in B z . When the scale length for such gradients is of order r, then this will move 
vertical magnetic flux radially at a speed comparable to the accretion velocity. This would 
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prevent radial infall from building up the vertical magnetic field at small radii to arbitrarily 
high values even in the absence of radial buoyancy. In gas pressure dominated disks it also 
prevents radial buoyancy from expelling vertical flux from the disk, even in the absence of 
a disk dynamo. In a radiation pressure dominated disk this effect will be smaller, and will 
allow a positive d r \nB z of order (P / P gas )r~ l . 

We conclude that the strong magnetic fields required for nonthermal emission from 
the central parts of AGN disks must be either be maintained by dynamo processes at 
small radii, or dragged in through a magnetically dominated, and relatively nonturbulent, 
outer disk. In the latter case, the vertical field will adjust to changes in the disk at a the 
characteristic rate of 

Tb 1 = £ ~ «f «. (31) 

The former case will allow the field to adjust on a dynamo time scale, which will be 
( IVishniac 1995b| ) 

T dynamo ~ ( 7T~) a ( 32 ) 

The nature of the dynamo at small radii is beyond the scope of this paper. Most 
work on disk dynamos has concentrated on the generation of toroidal fields, but the same 
processes will produce a poloidal field, albeit at a much reduced rate. In this case the field 
lines will tend to open up as differential rotation pumps energy into the disk magnetosphere, 
ultimately producing a large scale poloidal field. 



3. Implications for the Variability of AGN 



We have examined the radial buoyancy of poloidal field lines embedded in an accretion 
disk. In general the field lines will move outward through buoyant forces at a rate 
comparable to the inward drift caused by accretion. This leaves the question of their 



-16- 



global evolution unsettled. However, in radiation pressure dominated disks radial buoyancy 
dominates over accretion unless the magnetic field is strong enough to suppress turbulence. 
It follows that the global magnetic field embedded in accretion disks around AGN is likely to 
have its origin in dynamo processes taking place at small radii. In reaching this conclusion 
we have assumed that the angular momentum transport in these disks is mediated by 
magnetic field instabilities. Since the enhanced vertical buoyancy of magnetic field lines 
in such disks will tend to decrease the efficiency of this process ( |Vishniac 1995 bp it seems 



possible that this assumption may be violated in the presence of hydrodynamic transport 
processes. As an example, we cite the possible role of internal waves, which normally give 
a ~ (H/r) 2 ( [Vishniac fc Diamond 19"89|) , a dimensionless viscosity which is unlikely to 
dominate in disks where the radiation pressure is negligible. 

What does this imply about the variability of nonthermal radiation from AGNs? If 
the poloidal field is generated from within the disk, and moves outward through turbulent 
reconnection and radial buoyancy, then both its generation rate and its outward migration 
rate will depend on conditions within the disk. Consequently, variations in the state of 
the disk will lead to changes in the basic parameters of the magnetosphere, which will 
lead directly to variations in its nonthermal radiation. Variations in the mass accretion 
rate in the inner disk will arise from instabilities operating at larger radii and will reflect 
the thermal rates at those radii, roughly aQ(r). Since the time scales for disk variations 
will approach the light travel time across the magnetosphere for the inner regions of disks 
accreting near the Eddington limit it will be necessary to include time derivatives in the 
formulae for the structure of the magnetosphere. 

As long as the poloidal field has only a weak effect on the mass flow within the disk, 
compared to internal processes driving angular momentum flow, the effects discussed here 
are incapable of driving disk instabilities. However, if the radiation pressure dominated 
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part of the disk is subject to local instabilities over its entire radius, then such instabilities 
will lead directly to variations in the structure of the magnetosphere and the nonthermal 
emission it drives. For example, thermal transition fronts will produce narrow annuli with 
a large local variations in d r M, driving the magnetic flux in the direction of the local mass 
flow. Since the disk instabilities will reflect the thermal and mass transport time scales 
for a wide range of radii, from the inner edge of the disk out to the radius where the 
radiation pressure of the disk loses its dominance, these time scales will also be present in 
the nonthermal emission. The upper limit to these time scales will depend on the details 
of the disk structure, but will clearly be much longer than orbital periods in the inner part 
of the disk. However, the shortest time scales for variations driven in this way will be 
comparable to the period given in Eq. ([!]). 

This work has been supported in part by NASA grants NAGW-2418 and NAG5-2773 
(ETV) and by the Basic Research Project 95-5200-002 of the Ministry of Science and 
Technology, Korea (SJP). 
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